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1 Introduction 

Studies on the formulation of physical theories on the Moyal plane were initiated in recent 
times by Doplicher et al. fl]. Interest in such algebras was also stimulated by the work of 
string theorists who encountered them in a certain decoupling limit PJ. 

The d-dimensional Groenewold-Moyal spacetime is an algebra generated by 

elements x M (// 6 [0,1,2, • - - , d — 1]) with the commutation relation 

\_Xfj ,, xf\ iO^l , (1) 

0 jiu being real constants antisymmetric in its indices. In the limit 0^ v = 0, Xq and xi 
are time- and space- coordinate functions. If x = (xq ,x) is a point of when 6 IIU = 0, 
then 

x 0 (x) = x 0 , x.i(x) = Xi. (2) 

Thus Xq, Xi are operators in Ae(IP d ) which become time and space coordinate functions 
when Of, v = 0. There is an extensive literature on the formulation of quantum held theories 
(qft's) on M.6»(IKd) and on their phenomenology 13]. The focus of much of this work is 
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on space-space noncommutativity (% 7 ^ 0, i,j G [1,2, ...d — 1]). But it is time-space 
noncommutativity (9q * 7 ^ 0 ) with its implications for causality and foundations of quantum 
theory which leads to strikingly new physics. The formulation of unitary qft's when @0i 7^ 0 
is nontrivial and was carefully done already by Doplicher et al. [I]. 

I 11 recent papers mu, Balachandran et al. formulated unitary quantum mechanics 
when 9 0i 7 ^ 0 basing themselves on the ideas of Doplicher et al. Consequences of spacetime 
noncommutativity in the quantum mechanics of atoms and molecules have been explored 
by Balachandran and Pinzul JfT] • 

Previous work on ^4e(M a! ) was focused on the formulation of quantum theory. Effects 
of noncommutativity on classical waves and particles have largely remained untreated in 
particular when 9^ 7 ^ 0 (see however flj). In this paper we discuss “classical” waves on 
^ 4.0 (K. rf ), assuming time-space noncommutativity (6^- = 0, 9 0i 7^0). 

The approach adopted in Doplicher et al. |T| and subsequently in Balachandran et 
al. [ 3 ] to study space-time noncommutativity is different from the string theory motivated 
studies in the literature due to Gomis and Mehen [J\ and other authors | 8 ], which found that 
field theories on noncommutative spacetime are perturbatively nonunitary. As explained in 
detail in Balachandran et al. J2], in the former approach, the amount r of time translation 
is not “coordinate time”, the eigenvalue of xq. For 9 = 0, these two could be identified, 
while for 9 7 ^ 0, Xq is an operator not commuting with and cannot be interchanged with 
r. The separation of eigenvalues of £0 from the amount of time translation is the central 
reason for the unitarity of the theories as formulated in Doplicher et al. and Balachandran 
et al mm- This is analogous to the situation in quantum mechanics, where if p is the 
momentum operator, spatial translation by amount £ implemented by exp (i£p) is not the 
eigenvalue of the position operator x. 

I 11 the algebraic approach (which is mandatory if 9 )W 7 ^ 0), waves are elements of the 
spacetime algebra. That is the case also for the commutative space-time C'°(M d ). The 
act of observation, such as measurement of mean intensity over the time interval T, is 
represented by a state on this algebra. In the following section, we describe this approach, 
valid equally for commutative and noncommutative algebras. Subsequently it is applied to 
interference for a double slit experiment for the algebra ^(M^). I 11 cases where a double 
image of a star is formed by a cosmic string, it causes interference as well which is affected 
by 9 0i . This phenomenon is examined in the final section. 

Novel phenomena are observed in interference when 7 ^ 0. For instance if the time 
of observation T is too small, then, as indicated before, one sees constant intensity and 
no interference on the screen. Interference returns for larger times, but it is shifted and 
deformed as a function of (9 — , w : frequency of the wave). The familiar 

interference pattern is recovered only when ^ —> 0 . 

2 Classical Waves and Particles on Algebras 

The algebra 4.£i(lR d ) has generators x M with relations [x M , x v \ = 9^ u = —9 Ufl being real 

constants. We assume that % = 0, and orient 9q = (#01 , ■ • - 6 W-i) in some direction n. 
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Thus, for us, 


[xi , Xj] = 0 , [x 0 , Xi] = iOrii , 9 G R, n ■ n — 1. (3) 

We can set 9 > 0. This does not entail loss of generality since 9 flips in sign when 
Xi —> — Xi. So 9 0 = 9n, 9 > 0. 

For 6 = 0, is the algebra (7 0 (M d ) of functions on M0 We first outline the 

algebraic approach for 9 = 0. It generalizes easily to 9 0 0. 

i. Classical Theory on Commutative Algebra: 

Let us first examine waves. They are fields 0 on spacetime so that 0 G *4.o(® d )- It is 
enough to consider scalar waves. Then 0(x) for x = (x 0 ,x) e ^4 0 (^ d ) is tl ie amplitude of 
the wave -0(x) at x. It is the solution of a wave equation such as 

f d o ~ = °' ( 4 ) 

The algebra ^4o(^ d ) contains not just 0, but functions of 0 as well. It is reasonable to 
assume that a general element a G ^4o(^ d ) has the Fourier representation 

a= I d d ka(k)e iko *°e iU , (5) 

where x M are the coordinate functions: x^x) = 

We can measure many attributes of a wave. For example, we can measure its mean 
intensity over the time interval [xo — \ , Xo + f-]. It is 

1 r x o+| 

/ = - / dx o | 0 (x o ,x)\ 2 . ( 6 ) 

T Jx o-f 

We want to interpret this measurement as the application of a state on a particular element 
of the algebra since states are defined also for noncommutative algebras. 

A state u on a ^-algebra A with unity 1 is a linear map GDI, 

—> C (7) 

which is positive 

u(a*a) > 0 ( 8 ) 

and normalized : 

w(l) = l. (9) 

Thus states define probabilities and u(a) is the mean value of a. 

Coming back to ©, for intensity, we associate the observable / e ^4o(^ d ) where 

J=| 0 | 2 . ( 10 ) 
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Measurement of the mean value of a G Mo(ffi d ) at x in the time-interval [xo — ^ , xq + 
is represented by the state u where 


uj(a) 


1 

T 



dx o d(xo , x). 


So u ; depends on xq, T and x. Then 


(ii) 


I = U0). (12) 

Thus to an observable, we assign an element a G Mo(IR d ) and to a measurement, a state uj 
on Mo(K d )- The result of this measurement of a is uo(a). 

A classical particle too can be described by a similar formalism. Instead of working 
with it is best to include momenta also and work with M 2d_1 . A point of M 2 ^ -1 is 
(xo,x,p) where p denotes momenta. The algebra is then Mo(M 2<i_1 ). If a G Mo(M 2d_1 ), 
a(xo,x,p) is the value of the observable a at time x 0 for a particle with position x and 
momentum p. For energy E G ^(IR 2 ^ 1 ) in a possibly time-dependent potential, we can 
have E(x 0 ,x,p) = ^ + V{x Q ,x). 

States too can be introduced. For example, define u by 

1 r x o+| 

u(a) = — dx 0 a(x 0 ,x ,p). (13) 

Jxo-^ 

So uj depends on xo, T, x and p. uj{E) is the mean energy in the time interval [xo~ ^ , Xo + f] 
for a particle at x with momentum p. 

We, however, will not pursue point-particle theory any further. 

ii. Classical Waves on Noncommutative Algebra Mo(lR d ) 

As remarked already, states can be defined also on Thus to carry the discussion 

forward, we must identify waves in say by wave equations, associate observables to 

waves and define suitable states. We will do so in the context of interference and diffraction 
for d < 4 in what follows. But we must emphasize one strikingly new feature of 9 ^ 0. Then 
since Xq and x* do not commute, by the uncertainty principle, we can not simultaneously 
localize time in an interval T and sharply localise spatial coordinates. So a state like uj in 
m with exactly the same features does not exist for 6 ^ 0. We can at best approximate 
it. 

We consider free massless scalar fields ifj G for d = 2,3,4. Such massless scalar 

fields obey the standard wave equation 

(<% - V 2 )^ = 0 (14) 

for 9 = 0. We must find its analogue for 6 ^ 0. 

For simplicity, we choose n = (1, 0, 0), if necessary by applying a spatial rotation on x,;. 
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Let 



Po = adxi 



ad a f3 := [a , f3\. 

(15) 

Then 

P 0 x 0 = i , PoXi = 0 , i > 1. 

(16) 

So Po substitutes for and 

we can identify do with — iP 0 : 



§° 

1 

1 

. 

(17) 

Similarly 

d\ — > iPi = --adx o 

U 

(18) 

while 

d a — >d a =iP a , a = 2,3, 

(19) 

d a in (19) being the conventional differentiations. So the 

noncommutative elementary 

massless wave equation is 

(PS - P; - PS)i> = o. 

( 20 ) 

It has plane wave solutions 


( 21 ) 

with a standard dispersion relation: 



w 2 — k 2 = 0 . 

( 22 ) 


The general solution is a superposition of plane waves. 

We note that for electromagnetic (EM) waves, (I2UI) receives corrections in powers of 
9, since in noncommutative spacetimes, the EM Lagrangian gives nonlinear equations of 
motion. Interestingly enough monochromatic plane waves of the form m are solutions to 
the nonlinear equations of motion to every order in 9 HD- But their superposition is not. 
Nevertheless, as the inclusion of this effect will give only higher order corrections in 9 to 
our results, they are not treated in this paper. Similarly the use of “covariant coordinates” 
(cf. and references therein) will not affect leading order results in 9 and hence we work 
with the standard noncommutative coordinates. 

a) d = 2 : 

The problem we examine is the interference of two plane waves with the same frequency. 
It can be generalized, but several essential points are well illustrated by this example. Thus 
we consider 


'ip = i> k + i >- k , 

j) ±k = e ±ik ^ e ~ ikS:o , k > 0 . 


( 23 ) 
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We see that the intensity 

7=H 2 (24) 

is an operator in „4.<j(M 2 ) . 

It is not possible to achieve a state with a sharp localization in position and time. 
Instead, we look for a state uj with a reasonable spatial localization around a point x\, it 
will be rather delocalized in time. We define oo = oj 1 in terms of a density matrix according 
to 

Tryd 


Uj{a) = 


Try(l) ’ 


where 


7 = ^t(x o - x 0 )5(£i - £i) i> T (x 0 - x 0 ) 
Here by 5 and ^t we mean the following operators: 


S(x i 

4>r{xo - x 0 ) = 


Xl ) = — dke ik ^~ Xl \ 


t T 
2 


T 
' 2 


dX <5(x 0 — xq — A) 


Let x \) be the eigenstate of x\ for the eigenvalue x\: 

x\\x']) = a/Jx)) . 

Then 

<5(ah — xi)|a^) = 8{x' 1 — xi)|x , 1 ), 

where on the right hand side stands an ordinary delta function. 
ipT is defined on eigenstates \x' {] ) of x 0 ; 

xo\x 0 ) — Xg|xo) , 

^ T (xo - Zo)K) = ^t(x o - x 0 )K) . 

Here ipr(0 is the characteristic function on the interval [— 7 rf"] 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


tMO = 


1 for |f| < | 
0 for |£| > 0 . 


(31) 


(32) 


It is easy to show that 7 is a positive operator and that is a state. Now \j)\ = 
Hence for 9 = 0, we can write 7 = S(x 1 — Xi)^t(xq ~ Xo). The corresponding w 7 describes 
an experiment at spatial location x\ which is averaged uniformly over the time interval 
[x 0 — , x 0 + -|]. For 6 ^ 0, 7 is an approximation to such an experiment. 

But for 9 0, 7 does not have sharp spatial localization. If it did, then for 

a = <5(xi - yi ), (33) 

we should get 

uj(a) = 5(x 1 - 2 / 1 ) . (34) 
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Instead we find (cf. Appendices 1 and 2 ) 


2 e 

u; 7 (aj = — 


7rT(xi- yi) 2 


sin" 


T(x i - yi) 
26 


(35) 


As | -mx), it approaches to <5(xi — yd) as it should. 

It is important to emphasize that in order to interpret the action of the state cu 7 on \iJj\ 2 
as the measurement of intensity at a given spatial location, say x\, it is necessary to be able 
to localize x± with a reasonable precision which becomes sharply localized as 6 —> 0. ay. 
with density matrix 7 does this job perfectly: It approximately localizes the point where 
the measurement is taken (c.f. equation (ESI)) and it has the correct commutative limit. 

A convenient way to calculate the traces such as those in ESI) is to use the coherent 
states. Let 

Xq A ix\ x Xq — ict\ r x 

,— , a' = -^=—, ~ 

V29 V29 

and 

1 \ — h=g(zaJ —za) 1 n\ 

| z) = ev® 1 ; |0), 

where a\z) — -^=|^) and {z\z) = 1. We have 


a = 


[a, a 1 } = 1 , 


(36) 

(37) 


(z\x^\z) = X^ , 


(38) 


where z = x 0 + ix\. 

We can now compute ESI) using the resolution of identity 


d 2 z 


47 t0 




(39) 


In particular we find that 


Try = 


d 2 '< 


AttO 


WW) = 


T 

2 Ve 


(40) 


Appendix 1 contains the details. 

Note that PHll diverges as 6 —> 0. That is because h(£i — x \) is not normalizable in the 
commutative limit just like a plane wave. 

The object / of our interest is the intensity as measured by cuy. 


We find, using coherent states for example, that 

Try\^>\ 2 _ J 2 [l + (l — cos2rc(a;i — 9w)] for 2 6w <T . 


Try 

as Appendix 3 shows. 


for 2 Ow > T 


(41) 


(42) 
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Figure 1 : Variation of intensity / as a function of wx i, for fixed ^ The plots are for 
r — T anc l Q w 2 = Clearly, as the ratio ^ gets smaller, it converges to the 

commutative result. 



Figure 2: Variation of intensity / as a function of y — 6. ft is plotted at x\ — 0, u — 
10" 7 m _1 for the values of T = 2.5 x 10 7 , 4 x 10 7 ,6.5 x 10 7 ?tz. Note that for 6 = 0, I = 4 
and for any given T, I takes the value 2 and becomes independent of 9 at 6 = ^ x 10' and 
thereon. 









This result is remarkable. It asserts that there is no interference at all if ^ > |/ Thus 
the higher the frequency, the larger is the time of observation needed to perceive interfer¬ 
ence. There is interference for ^ but its pattern is shifted and extrema modified 

depending on frequency and time of observation. Figure 1 illustrates the phenomenon. 
We recover the usual pattern when ^ > 0, and in particular in the commutative limit. 

Variation of / w.r.t. 2 9w is plotted in Figure 2. 

b) d = 3: 

We now consider the algebra ^4e(M 3 ). A good problem to study is Young’s double slit 
experiment. For this, let us imagine a screen (a line) along the 2-direction. The two slits 
are a distance a apart and the line joining them is also parallel to 2-axis. This line is at a 
fixed distance from the screen. 

As spatial coordinates commute, there is no intrinsic difficulty in realizing an arrange¬ 
ment like the above with arbitrary precision. 

Let n= (1,0). Then time-space non-commutativity is expressed as 

[x 0 , £i] = id , [x 0 , x 2 \ = 0 . (43) 

The wave is considered to be associated with a massless scalar field. Then as in (23), 
a plane wave is 

■01 = e ik ' r e ~ lkx °, ? = (xi, x 2 ), k • ? = k t Xi , |k| = k — \Jkj + k% . (44) 

We can also spatially translate this wave by a and evolve it for time r by applying 
e ),p-a e —* p 0 t, the result is still a plane wave: 

e -iP 0 T e iP-a r e ~ikx 0 ^ _ ■?g-ifczo) g* ka g-* fer j'qgj 

The last factor is the complex-valued phase of the wave. It lets us to unambiguously 
compare the phases of waves related by space-time translations. As the phase does not 
change if a — |r for k 0 , the phase being zero for k = 0 , we can say that the wave 
travels in the direction k as in the 9 = 0 limit. 

We want to consider the interference of the waves from the two slits at a point P on 
the screen, assuming for simplicity that they have the same frequency w. 

Let k and k 7 be the directions of propagation from the slits to P. Then the wave at P 
is 


+ 4>2 , 

= e * kf e ~* fc£ o 5 = e* k '- t' e -ikx 0 ; 


(46) 


where 

k = kik , k ; = k'\t , r' = r + a (47) 

Here a = a(0,1) is the displacement of the primed slit relative to the other one. 
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Note that the waves do not acquire phases as they arrive at P from the slits a time r 
later as the remark after © shows. 

We generalize the density matrix 7 according to 

f = tp T (xo - £o)f>(xi - xi)5(x 2 - x 2 )i> T {x 0 - x 0 ), (48) 


where <5(ay — ay) and are given by (1271) and (I2H1) respectively, while S(x 2 — x 2 ) is a 
regularized version of the delta function centered at .x 2 (see below). Consequently, the state 
generalizing © is given by 


(Up(d) 


TrYa 

Trf 


(49) 


Then the intensity / due to -0 at P is given by cap ('i/4?/;). 

Note that the standard 5-function 5(x 2 — x 2 ) is not normalizable, having infinite trace. 

Hence, to regularize its contribution to the traces, we replace it, for example by S(x 2 —x 2 ) = 
|a)(a| where 

(x 2 \a)(a\x 2 ) (50) 


is a Gaussian of width d centred at x' 2 = x 2 . For this purpose we introduce the momentum 
P 2 conjugate to x 2 with the eigenfunctions e ip '^ x ' 2 ~ X2 ^ and eigenvalues p' 2 € (— 00 , 00 ). 
Consider 


{x' 2 \ a) = 


1 yfd 


V^F^/4 


dp 2 


Jp' 2 {x' 2 -X2) -±d 2 p% 


7T l / A \fd 




(51) 


m will do the job: © is a Gaussian, which in the d —» 0 limit is a delta function centred 
at x 2 , while Tr|a)(o;| = 1. 

With 6(x 2 — x 2 ) = |a)(a|, traces can easily be computed in the basis | z') <g) \x 2 ) (| z') 
being the coherent states used in the previous section). We observe that TrY = while 
for the intensity / we hnd 


/ = u= lim 2 

d^O 


1+1 


^— M) e -id 2 (k2-k^ cos (_ ^)( Xl - k9) 


= 2 


1+1 


9\ki — k\ 

¥ 


—{k 2 - k 2 )x 2 - k 2 a 

1 ^ cos ^(Ay — Ay) (a; 1 — k6) + ( k 2 — k' 2 )x 2 — k 2 a 

for 6 1 Ay — Ay \ < T , 


1 = 2 for 8\ki - k[\ >T. 
Note that the final result is independent of d. 


(52) 
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0 


identify 



Figure 3: Double image of Q is observed by 0 = 0'. 


This result is similar to the d =2 case. There is no interference at all for 9\ki — k\ > T 
and the interference pattern is distorted for 9\ki — k[\ < T in a similar manner encountered 
in (gH) . 

When 0\k\ — k\ < T, the dependence of / on ( k\ — k[)x± can be plotted for fixed values 
of 6 ^ kl T — . With suitable choice of the phases in K52D, the plots will look similar to those 
in Figure 1. 

c) d =4 : Interference Phenomena from Cosmic Strings 

As an example we here study interference of waves from a distant source caused by a 
cosmic string. For simplicity we consider a straight cosmic string. The metric around it is 
flat with deficit angle 87rG^t, where /i is the mass per unit length of the string and G is the 
gr avit at ional const ant. 

Figure 3 shows a spatial slice where Q is the source and S is the string. We assume that 
S is normal to the spatial slice as it is sufficient for our purposes. The string causes deficit 
angles and requires the straight lines SO and SO' to be identified. This identification 
causes a double image of Q. A = OSO’ is the deficit angle 8nG/a. Q, S , O and O' are on 
a plane V. Calling l the distance from the string to the source and d the distance from the 
observer to the string, the angular separation of the images is |]2| 

E = £ + if = (53) 

(where £ and r/ are explained by Figure 3.) 

We note that all the spatial coordinates commute with each other. Hence for the set-up 
above there is no intrinsic difficulty in defining the spatial directions and their orthogonality. 
They are exactly the same as in the commutative case. 

ft should be clear that the problem is effectively of two spatial dimensions. On the 
spatial slice any fixed vector from the string can be taken as the direction noncommuting 
with the time coordiante. 
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Figure 4: The observer e at 0 = 0' can barely observe a double image of the source Q. 
Keeping her distance from the string fixed, she can shift her position along the arc A and 

can observe a double image of Q from O = O'. At 0 = 0' she again can only barely 
observe a double image of Q. 


Consider now the identified observers at O and O'. Call her e. It is not hard to see 
from Figure 4 that keeping her distance d from the string S fixed, she can observe double 
images due to Q on all points on an arc A of length dA. (See Figure 4). 

Suppose now that she can observe the intensity of waves of definite frequency w coming 
from Q. Let us denote the wave vectors of the two plane waves and ^2 emerging from 
Q that reach the observer e at O = O' by k and k'. respectively. At 0 = 0' the angle 
between k and k' is equal to S. The intensity observed by e at O = O' is then given 
by (JS2J), where a in that formula (the separation between the slits) is now zero as there 
is only a single source Q. The observer e can shift the observed intensity by shifting her 
position on the arc A. Thus we can think of A as a screen, where the interference pattern 
is recorded. 

After some trigonometry we find that (see Appendix 4) 


| h - k[ 


{k 2 + k' 2 ) tan 



(54) 


From this result it is possible to estimate an upper bound on |) in order that the two 
waves interfere. Recall that for this we must have e ^ kl T < 1. Substituting for \k± — k[\ 
from we get ^AbA < | co t = |. It is known that double images of quasars are usually 
separated by a few arc seconds [Tl] . Taking for example an angular separation of 5 arc 
seconds puts the bound e \ k2 + k ^ < 8.25 x 10 4 . We note that k 2 + k' 2 = ak where 0 < a < 2. 
Therefore, ^ < 8,2 ^, 1C|4 and for a given k the lowest upper bound will be approached as 
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a —> 2. Finally, suppose that a wavelength at the red end of the visible spectrum is 
observed, say at A = ^ = 700?rm. Then we find ^ < 4.6 x 10~ 3 m. Only under this 
condition is the interference observable for light with wavelength A = 700nm. 

3 Conclusions 

In this work we have studied the general theory of waves in Groenewold-Moyal spacetimes 
where time and space coordinates do not commute. We have given the rules for the 
measurement of their intensity and applied them to interference and diffraction phenomena 
in spacetimes of dimensions d < 4. The latter produced novel physical results. Namely, we 
found out that for observation times T which are so brief that T < 29w, no interference 
can be observed. For larger times, the interference pattern is deformed and depends on 
It approaches the commutative pattern only when ^ —> 0. These results are given 
concretely by the equations <H2D and (|52p for d = 2 and d = 3, respectively. Finally, we 
have used these results to discuss the interference of stellar light due to cosmic strings, 
where with the help of the present stellar data we have estimated that for a given k, we 
must have ^ < 4x40 to observe interference. 
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Appendix 1 

In this Appendix we calculate Try where 7 is given by (26). In the coherent state basis 
| z') we have 
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Using (27) and (28) in flA. II) . we get 


Try = J dX\dX 2 Tr 6(x 0 — xq — Ai) h(fi — xi)5(fo — x o ~ A 2 ) 

' ' dX^dhdhdh Tr e ifci(*o-*o-Ai) e ife(ii-* 1 ) e ifc 3 ( 4 o-*o-A 3 ) (A.2) 

dAidA 2 dfcidA: 2 dA:3e i[0fe2fe3 - fel(;ro+Al) - fe2X1 - fc3{xo+A2)1 Tre i(fel+fe3)£o e ife2£l 


(2tt) 3 

1 

(2^)3 


Using (36) and in the coherent state basis we have 

2y e *( fe i+ fe 3)£ogifc2ii _ / z /| e i(fe i +fc 3 )io e ifc2xii r /\ 


e— x ! 


_ f dx odx' 1 ,e |(fc 1+ fc 3 )2 + fc2}-i» ( fc 1+ fc 3 )fc 2 life+fcsK ifc 2 

./ 2 tt 0 

= ^ e [!{( fc l+ fc 3) 2 +fcl}-i|(fcl+fe3)fe2]^^ i + ^ 

C? 

= + fc 3 ) 6(fo), 


(A.3) 


where we have used z' — x' 0 + ix\ and d 2 z' = 2dx , 0 dx , 1 . Substituting (I A. 3D in (A.2) and 
remembering that — ? < A * < ?, z = 1,2, we get 


Tr 7 = 


4n 2 0 

1 


2tt9 

T 

2^6 


I dA 1 dA 2 d/c 3 e* fc3(Al+A2 ' ) 
dAidA 2 h(Ai + A 2 ) 


(A.4) 


Appendix 2 


I 11 this Appendix we give the derivation of (35). Using (26), (27), (28) and (33) we get 
Tryd = Tri) T {x 0 - x 0 )S(xi - x^t^o - x 0 )5(x 1 - 


1 


(2tt)4 

1 


(27r) 4 
idA 2 

coherent state basis we get 


d/l 7 Ye* fcl Uo-zo- A l) e *fc 2 (zi- 3 : l) e * fc lUo-a:o- A 2 ) e ifc 2 (£l- 2 /l) m 5 ^ 

d/J, ^ 9k2k3 ~ k d x 0+M)-k2Xi-k 3 (x 0 +\2)-k4,y] rp^ e i(k 1 +k 3 )x 0 e i(k2+k A )x\ 


where dfi = dX\dX 2 dk\dk 2 dkzdk± and —7 < A * < 7 , i = 1 , 2 . Using (36) and in the 


r p re i{k 1 +k 3 )x 0 e i{k 2 +k 4 )xi 


^ e |[(fcl +fc 3 ) 2 +(fe2+fe4) 2 ] 5 ( A;i + jfc 3 ) ( J(fc 2 + fc 4 ) 

0 

2 tt 

—h(fci + k 3 )5(k 2 + h) (A. 6 ) 
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where we have used a technique similar to that was used to derive (IA.3f) . Using (IA. 6 I) in 
(A.5) we get 


Trja = 




d\id\ 2 e 1 9 e 2 9 


tt 2 (xi - ?/i) 2 

Using (61) and (58) in (25) we get 

2 1 


Sill' 


= 


ttT (xi - 2 /i) : 


Sill 


^(^i - yi) 

29 


Tjx i ~ yi) 
29 


(A.7) 


(A. 8 ) 


Appendix 3 

In this Appendix we provide the derivation of the intensity formula (42). Using (23), 
we get 

Try|V)| 2 = 27>7 + e ~^ Tr ^k^ + e ^ Tr ^ e 


= 2 


Try + Re(e _2ifc e Tr^e 2ikil ) 


(A.9) 


Substituting for 7 from (26), we can write 


Trje 


2ikx\ _ 


(2 n)‘ 


d^ ie -i[k2k39+ki( x o+\i)+k2 x i+k3( x o+\2)\ r p r ^i{ki+k3)xo^i{k2+2k)x\ 


where d/r = dXidX 2 dkidk 2 dk 3 . Using (36) and in the coherent state basis we get 

Tre i(k 1 +k 3 )xo e i(k 2 +2k)x 1 = ^ e lp 1 +fe 3 ) 2 -(fc 2 + 2 fc) 2 ]^( fcl + k 3 )6(k 2 + 2k) 

0 


2n 

T 


5{k\ + k 3 )5{k 2 + 2k). 


(A.ll) 


Hence 


p2i(kx\ — k?0) 

-2ikH TrAfe 2ikxx = / d\ 2 / dAi <5(Ai - A 2 - 29k) 


(2it)9 


O 2i(kx\—k 2 0) C~ 2~2 kO 


( 2 . 71)9 


d\ 


T 
’ 2 


2 - 


(A.12) 


Remembering that 9k > 0, we see that the integrand in (1A.12I) vanishes if 29k > T. If 
29k < T instead, we have 


p 2i(kxi—k 2 6) /*?—2 kd 
^—Zik^Orp^z^likxi _ J \ 

Trie ~ ( 2 Tje J_r ix - 


^2i{kx\ — k 2 6) / rp 

7 t9 


-9k 


(A.13) 
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0 = 0 ' 


Figure A. 1: The diagram provides the definition of angles ifj, i(' and (3. 


From (1A.9[) . (IA.12H . (I A. 131) and (I A. 411 we finally get 

Trj\'ip\ 2 _ f 2 [l + (l — cos2zr>(xi — 9w)] for 2 6w < T, 
' ' 2 for 2 9w>T 


I = 

Tr'y 

where we have used w = k. 


(A.14) 


Appendix 4 


In this appendix we give a derivation of the result PI). From the figure above we see 
the following relations between the angles 

ijj' — -0 = H, (A.15) 

2/3 + E = 180 =y (3 = 90 - |, (A.16) 

2/3 + ip' - ip = 180 => i(j' + (3= 180-/3 + ^- (A.17) 


Noting that cos^ = jf , sin ip = , 


n.(kxr) 
kr 

k x r) and \k\ = k' = k, we explicitly have 

1 


(n being the unit normal in the direction of 


cos^=- (h- + k 2 —) sin-0 = - (h— - k 2 -) , (A.18) 

k V r r / k \ r r / 


and similarly 


cos?//' = 


*4 


Xi 




(A.19) 
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Using the above relations we find that 


cos (/? + ip') = 


k 


, . . _ . _ > X]_ I ... ZL . £1 \ Xi 

A:-, sin —hfc, cos —-h re 9 sin- tc, cos — — 

1 2 2 2 Jr V 2 2 1 2 Jr 


cos(180 — (3 + ip) — — 
k 


( —k\ sin A 

+ k -2 cos A 

— + 1 

k >2 sin A 

— k\ cos — 

X2 

Lv 2 

2^ 


V 2 

2 ) 



(A. 20 ) 


Due to (1A.17I) the difference of the the cosines above vanish. For this to happen, we observe 
that the coefficients of — and — must separately vanish. Thus the coefficients of — must 

J' rp A V Y" 

satisfy 

(k 2 + k' 2 ) sin j = (k[ - ki) cos A , (A. 21 ) 

and hence 0 . 
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